SPHERES, Jiilich's High-Flux Neutron Backscattering Spectrometer at FRM II 
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I. INTRODUCTION 

Neutron backscattering is a versatile technique for 
measuring the dynamics of spins, atoms, and molecules 
on a GHz scale. Typical applications include hyperfine 
interactions in magnetic materials, molecular rotations, 
diffusion, and relaxation in complex systems. ^'^ 

Historically, three instrument generations can be dis- 
tinguished. Proposed by Maier-Leibnitz in 1966,^ neu- 
tron backscattering was first realized in a test setup at the 
4 MW reactor FRM (Garching, Germany). Count rates 
were of the order of one per hour, and the signal-to- noise 
ratio was 2:1 at best so that it took several weeks to mea- 
sure a single spectrum^ Following this demonstration of 
principle, first-generation backscattering spectrometers 
were build at Jiilich ( "7r-Spektrometer"— ) and Grenoble 
(INlOa-8, IN13^); they were reviewed in Refs.[llMl 

In the second- generation instrument IN 16 of Institut 
Laue-Langevin (ILL) at Grenoble, consequent use has 
been made of focussing optics; a deflector chopper allows 
to operate a spherical monochromator in exact backscat- 
tering geometryfi^^— Another instrument of this type is 
planned to be built at the Australian research reactor 
OPALJ^ 

In third-generation spectrometers, the focussing op- 
tics is made even more efficient by a phase-space trans- 
form (PST) chopper. This device, invented by Schelten 
and Alefeld in 1984fil enhances the neutron spectrum 
available at the monochromator, and allows for a par- 
ticularly compact instrument design. The first realisa- 
tion has been the High-Flux Backscattering Spectrome- 
ter (HFBS) at the National Institute of Standards and 
Technology (NIST, Gaithersburg, Md.).i8.i9 

In this article, we present SPHERES, the SPcctro- 
meter for High Energy RESolution built and oper- 
ated by JCNS (Jiilich Centre for Neutron Science) at 
FRM II (Forschungsneutronenquelle Heinz Maier-Leib- 
nitz of Technische Universitat Miinchen at Garching, 
Germany). A third PST based backscattering instru- 
ment is currently under construction at the ILL (project 
IN16B) where it will replace both INIO and IN1632i 

The development of SPHERES has been documented 



in several publicationsi^^— First productive experiments 
were performed in June 2007, and routine operation 
started in April 2008. Since then, two thirds of the beam 
time is offered to external users. Experiment proposals 
are invited twice a year through the JCNS user office. 
By now (March 2012), nearly one hundred experiments 
have been performed, and more than twenty publications 
have appeared. In the following, we describe the design 
of SPHERES, its hardware and software components, 
and its performance characteristics. The supplemental 
material contains a list of suppliers and some technical 
details.^ 



II. INSTRUMENT DESIGN 

A. Instrument Layout 

In a backscattering spectrometer, strong energy filter- 
ing of incident and scattered neutrons is achieved by 
Bragg reflection from perfect monochromator and ana- 
lyzer crystals under angles close to 180°. This peculiar 
geometry leads to some complications in the spectrome- 
ter design; it entails the need for a primary beam deflec- 
tor and a duty-cycle chopper. A particularly compact 
design is achieved if these two functions are realized by 
one and the same component: by a chopper that bears 
deflector crystals on its circumference. As an important 
side-effect, the Bragg deflection by moving crystals re- 
sults in a PST that enhances the useable flux at the 
monochromator. 

The resulting instrument layout is shown in Fig.[T] It is 
similar to the layout of HFBSr^ with three major differ- 
ences: (1) HFBS is operated under vacuum whereas the 
SPHERES housing is gas tight but not vacuum proof; 
(2) at HFBS the monochromator is located outside the 
instrument housing; (3) in the secondary spectrometer 
detectors and analyzers have changed sides with respect 
to the primary beam, resulting in a larger scattering an- 
gle range at SPHERES. 
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FIG. 1: Basic layout of SPHERES. 



Bragg reflection for energy selection 
final neutron wavelength Af (A) 
final neutron energy E{ (meV) 
final neutron velocity vt (m/s) 
Bragg reflection for deflection 
deflection angle 2qo (°) 
Schelten-Alefeld velocity vq (m/s) 



Si(flf) Si(3fl) 

6.271 3.275 

2.080 7.626 

630.9 1208 

PG(002) PG(004) 

41.65 25.14 

300.1 290.5 



TABLE I: Parameters that depend on the choice of Bragg re- 
flections for energy selection and for deflection. The current 
configuration of SPHERES uses silicon (111) in monochroma- 
tor and analyzers and pyrolitic graphite (002) as deflector. 



B. Crystal Choices 

The basic design decision for a backscattering spectro- 
meter is the choice of the monochromator and analyzer 
crystals. Initially, SPHERES was foreseen to be con- 
vertible between a Si(lll) and a Si(311) configuration. 
Si(lll) is the standard at most backscattering spectro- 
meters. The shorter wavelength of Si(311) gives access 
to a larger range in the scattering wave number q, with 
promising application niches in the study of diffusion on 
very short length scales.— However, the lower reso- 
lution (which goes with A~^) and background problems 
from Bragg reflections (for instance from sample holders 
or from cryostat walls) have so far prevented wide-spread 
useage. In commissioning SPHERES, the original me- 
chanics had to be modified in many details in order to 
block neutronic background channels and to provide bio- 
logical shielding. As a result, the instrument has lost its 
convertibility; the Si(311) option has been abandoned for 
the foreseeable future. 

Another important choice is pyrolitic graphite (PG) for 
the deflector crystals. Once deflector, monochromator, 
and analzyer crystals are chosen, the neutronic parame- 
ters listed in Table U are fixed. The angle 2ao between 
incoming and deflected beam is obtained from the Bragg 
equation 2cosq;o = ^/{c/l) where c = 6.710 A is the 
lattice constant of PG, and I the order of the Bragg re- 
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FIG. 2: Neutron path versus time-of-flight: (a) As realized 
in SPHERES, with chopper opening frequency /i = Vf/AR. 
(b) Originally planned for SPHERES (and realized in HFBS), 
with /2 = 3/i . The upward arrows describe neutrons that are 
directly scattered from the sample into the detector. 
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C. Time-of-Flight Logic 

Another decisive hardware choice is the radius R — 
2000 mm of the analyzer spheres. This parameter fixes 
the time of flight t :— 2R/vf — 6.34 ms for a round trip 
from the sample to the analyzer and back. During this 
round trip, the detectors see n pulses of directly scattered 
neutrons, interlaced with n — 1 pulses of energy- analyzed 
neutrons. If both types of pulses have the same duration 
T„, then r„ = r/(2n — 1), with n — 1,2,... Accordingly, 
the chopper must open with a frequency /„ = 1/(2t„). 

In the primary spectrometer, neutrons are deflected 
by the closed chopper. After a round trip to the mono- 
chromator and back, they must hit the open chopper 
phase. Therefore, the distance chopper - monochroma- 
tor LcM must obye 2Lcm/v{ — (2m — l)r„. Since the 
monochromator should have about the same curvature 
as the analyzers, the only practicable choice is Lcm — R 
and m = n. All this is illustrated by the time-distance 
diagrams in Fig. [2l 

HFBS uses n = 2, and so will IN16B. The present 
chopper of SPHERES also was originally designed for 
n — 2, but because of mechanical instabilities it is only 
operated with n = 1 . This reduces the efficiency of the 
PST, but it has a positive effect upon the signal-to-noise 
ratio, to be explained below (Sect. lyTCl) . 



III. INSTRUMENT COMPONENTS 
A. Neutron Guide and Velocity Selector 

SPHERES is located in the neutron guide hall West 
of FRM II. Like all instruments in this hall, it re- 
ceives neutrons through beam tube SRI from the cold 
source. SPHERES uses an end position at the sec- 
ondary neutron guide NL 6a. This guide has a cross- 
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FIG. 3: Deflector crystal assembly in the PST chopper. In 
both drawings, the neutrons come from the left. Left draw- 
ing: cross section. Right drawing: view from outside towards 
the chopper axis (the closing Al element and the screws are 
not shown). The deflector stacks consist of three pyrolitic 
graphite mosaic crystals, tilted by aluminum wedges, and fol- 
lowed by a ^''B4C absorber plate. 

section of 60 x 120 mm^; its curved segments have a 
radius of 1000 m. The guide is shared with two up- 
stream instruments, MIRA^S and DNS, that use thin 
crystal monochromators to extract neutrons at shorter 
wavelengths than those needed at SPHERES. 

The instrument shutter of SPHERES and its mechan- 
ical neutron velocity selector are located in two separate 
biological shieldings. The velocity selector (EADS As- 
trium) reduces the incoming neutron spectrum to about 
6.27±0.34 A. It is followed by a convergent neutron guide 
("anti-trumpet", by S-DH) that reduces the beam sec- 
tion by a factor of 11.5 to 25 x 25 mm^. Neutrons leave 
this guide system with a divergence of about ±47 mrad 
horizontal, ±64 mrad vertical. For details on all this, see 
the Supplemental Materialj^^ 



B. Phase-Space Transform Chopper 

The phase-space transform (PST) chopper is the key 
component of a third-generation backscattering spectro- 
meter. As said above, it combines the functionalities of 
a beam deflector and of a duty-cycle chopper. Half of 
its circumference bears deflector crystals that redirect in- 
coming neutrons towards the monochromator. The other 
half circumference is open. When the chopper is open, 
neutrons coming out of the neutron guide are transmit- 
ted towards a beam stop that is integrated in the chopper 
housing, and neutrons coming back from the monochro- 
mator are transmitted towards the sample. 

To enable spectral scans, a sufficiently large wave- 
length band must be forwarded towards the monochro- 
mator. This requires deflector crystals with large mosaic- 
ity. Pyrolitic graphite (PG) crystals, specified for a mo- 
saicity of 2.5° fwhm (Advanced Ceramics) were selected 



(a) 79 Hz, 300 m/s (b) 26 Hz, 1 00 m/s 




FIG. 4: Chopper operation mode, (a) originally intended, (b) 
currently realized. The white square shows the illuminated 
area. The arrows indicate the rotation during r = 6.34 ms; in 
this time, deflected neutrons travel from the chopper to the 
monochromator and back. 

according to rocking curves measured one-by-one on the 
triple-axis neutron spectrometer UNIDAS at Jiilich. To 
further enhance the horizontal mosaicity, the crystals are 
assembled in stacks of three, with tilts of 0°, +2.5°, and 
—2.5° imposed by aluminum wedges (Fig. 

Intensity and width of the forwarded neutron band are 
further enhanced by the horizontal motion of the defiec- 
tor crystals. This phase-space transform 'from white to 
wide' has first been suggested in a conference proceedings 
of limited circulationjii an easily accessible summary has 
been provided in Ref. [l^. The Bragg reflection by mov- 
ing mosaic crystals compresses the energy distribution 
and thereby enhances the flux in the acceptance range of 
the backscattering monochromator. In accordance with 
Liouville's theorem, the compression of the energy dis- 
tribution is accompanied by a widening of the angular 
distribution. This is not a severe drawback as long as 
the angular resolution of the entire instrument is limited 
by that of the secondary spectrometer. 

Assuming an almost coUimated incoming beam, an op- 
timum transform is predicted for a deflector speed of 

sinao , , 

^^0 = Vi — 1) 

cos 2ao 

with the neutron velocity Vf and the PG deflection an- 
gle ao as given in Tab. H For Si(lll) and PG(002), 
one obtains vq = 300.1 m/s. In the meantime, several 
simulationsiSi^i^^ and a measurement^^ indicate that un- 
der realistic conditions the optimum deflector speed is 
only about 250 m/s. 

The current chopper rotor has a radius of 650 mm; the 
crystals are located at a mean radius of i? = 600 mm. 
The rotor has A: = 3 crystal bearing wings of 60° each. 
Originally, the chopper was intended to operate with a 
pulse fragmentation index of n = 2 (Fig. 2^): Neutrons 
that have made the round trip chopper - monochromator 
- chopper would hit the second but next opening follow- 
ing the wing by which they had been deflected. The 
rotation frequency /lot = (2n — l)fn/k = f2 — 78.9 Hz 
would result in a deflector speed of 27ri?/i.ot = 297m/s. 
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Tensile forces for such R and /rot are considerable. 
To keep them within admissible limits, the high-strength 
aluminum alloy alumold (ALZNMg2Mn) was chosen, and 
the geometry of the disk was optimized with finite- 
element calculations, resulting in a construction with 
three spokes and two recesses for each of the three 
wings, 

Thermal instabilities of the bearings and mechanical 
resonances prevent us from running the current rotor 
with the intended 78.9 Hz. As a fallback solution, we 
are operating the chopper with n — 1 and /rot = /i/3 = 
26.3 Hz (Fig. |4)d). A new rotor is currently under con- 
struction; it will have an asymmetric rotor {k = 1) to 
achieve a higher crystal speed while maintaining n = 1. 
As an engineering compromise, a defiector velocity of 
225 m/s has been chosen, somewhat below the optimum, 
which, according to the simulations, is quite fiat. 



C. Monochromator and Analyzers 

In the standard configuration, the Si(lll) reflection is 
used in monochromator and analyzers. Single-crystalline 
hexagonal wafers with a side length of 58.5 mm and a 
thickness of 500 |J.m were glued onto spherical supports 
with a curvature radius of i? = 2000 mm. The bend- 
ing stress induces a lattice-constant gradient that signif- 
icantly improves the backscattered intensity for a mod- 
erate cost in resolution. 

The latter can be estimated quantitatively as follows. 
The relative change of the lattice constant at a distance d 
from the neutral surface is nd/R with a material constant 
K = 0.776 for Si(lll)i^ Assuming a rectangular distri- 
bution of lattice constants, and neglecting Darwin tails 
and other dynamical scattering effects, one would expect 
a rectangular distribution of backscattered energies with 
a full width of AE = 0.81 |J.eV. Convoluting this dis- 
tribution with itself, to account for monochromator and 
analyzer, one obtains a triangular distribution with the 
same full width AE. This is in reasonable agreement 
with the observed overall resolution of SPHERES, which 
is in first approximation a Gaussian with a full width at 
half maximum of about 0.62 |j,eV (Sect. IVID|) . 

An alternate monochromator/analyzer set uses pol- 
ished wafers with a cut-in grid supposed to make them 
stress free. The resolution is about 0.4 |xeV, but the fiux 
is so much lower that this configuration is of little use, 
except possibly for the study of narrow tunneling or hy- 
perfine splitting. 

Additional monochromators offer the possibility of 
measuring with a fixed energy offset. A CaF2(lll) mono- 
chromator, providing an offset of 24 |j,eV with respect to 
the Si(lll) analyzers, is already available. 

The monochromator support, made of carbon- fiber re- 
inforced plastics, has a size of 500x250 mm^ and a spher- 
ical curvature of about R + ics/2 = 2180 mm, where 
Leg = 360 mm is the distance between the chopper and 
the nominal sample location, so that the monochroma- 



tor acts as a hollow mirror that maps the chopper's PG 
crystals onto the sample. 

The small-angle analyzers consist of eight massive alu- 
minum rings mounted on a common support. Each ring 
has a width of 95 mm; between two rings there is a gap 
of 10 mm. The innermost ring, with an inner radius of 
165 mm, proved useless because it receives fringe inten- 
sity from the direct beam; therefore it was covered with 
Cd. The other seven rings cover scattering angles 2i? 
from 9.0 ± 1.4° to 27 ± 1.4°. Typically, two rings focus 
into one detector. 

The large-angle analyzers consist of seven shells, made 
of cast aluminum and finished by turning They have a 
height of 2000 mm, and a maximum width (in the scat- 
tering plane) of about 530 mm. One shell covers a solid 
angle of about 60° x 15° ~ 0.27 rad^. The seven shells 
are located at scattering angles from 30° to 135°, with 
only small gaps between them. There is no one-to-one 
correspondence between individual shells and detectors; 
the seven shells must rather be seen as one contiguous 
surface. 



D. Monochromator Doppler Drive 

A novel, resonance-free monochromator drive has been 
developed jointly by Forschungszentrum Jiilich and an in- 
dustrial company ( Aerolas) . The development costs were 
shared with the ILL, which received an identical drive for 
IN16 and for future use with IN16B. The monochroma- 
tor shell is mounted on an arbor that is supported by two 
air bearing collars. Driven by a linear motor, it performs 
a sinusoidal motion. The displacement amplitude Z can 
be chosen between ±25 and ±75 mm. The maximum 
frequency is about _F = 10 Hz, resulting in a maximum 
velocity amplitude oi V = 2ttFZ = 4.7 m/s. 

To give an idea of the engineering difficulties, let us 
note that the peak acceleration is 2'kFV — 295 m/s^. 
The moving mass is nearly 10 kg, and the maximum me- 
chanical power is nearly 7 kW. The linear motors gen- 
erate a heat fiux of up to 2.5 kW that is removed by 
water cooling. The peak current of each of the two mo- 
tors is about 70 A at 600 V; much of this electric power 
is recovered during deceleration. 

We define the monochromator velocity v to take pos- 
itive values if the monochromator is moving towards 
the sample. In the rest frame of the monochromator, 
backscattered neutrons have the velocity Vf = 630.9 m/s. 
In the laboratory frame, they have the energy 

E,{v)^^{v, + vf. (2) 

The analyzers select scattered neutrons with energy Ef = 
i?i(0). In consequence, the sample energy gain in the 
scattering process is 

huj = E;-Ef= (2- + -4) Ef. (3) 

\ f f Vf^ J 
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The quadratic term causes a small asymmetry in the dy- 
namic range. With V = 4.7 m/s, accessible energy trans- 
fers hco extend from —30.9 to 31.2 |j^eV. 

E. Backscattering Detectors 

The detector bank for large scattering angles consists 
of eleven counter tubes 24NH15, bought from Eurisys, 
filled with 4 bar of '^He. They have a diameter of 1 inch 
and a height of 210 mm. 

The tubes are mounted vertically. They are alternately 
arranged on two circles with radii of 110 and 130 mm 
at intervals 10°, so that complete coverage is achieved 
for an angular range of 110° in the scattering plane. In 
front of the tubes (starting at a circle with radius 80 mm 
around the ideal sample centre) there are radial, vertical 
Cd collimator plates. The vertical opening of the detector 
block has been reduced to about 80 mm in order to block 
neutrons that do not come from the analyzers. 

For small scattering angles, front-window counters 
10NH3/5X P (also from Eurisys) are used, filled with 
3.5 bar of ^He. Their diameter is 2 inch, the active length 
is 1 inch. They are mounted around the exit nose of the 
chopper housing. 

F. Diffraction Detectors 

SPHERES is equipped with position-sensitive detector 
tubes (Reuter-Stokes) for auxiliary diffraction measure- 
ments to be performed simultaneously with the spectral 
scans. The tubes are located in the horizontal scattering 
plane in front of the analyzers; inevitably, they shield a 
few percent of the effective analyzer surface. At present, 
six tubes, with a diameter of 1 inch and 535 mm of active 
length, provide coverage of scattering angles from about 
30 to 120°. 

The resolution of this diffractometer is limited by the 
horizontal spread of incident neutron directions, which 
in the worst case is ±(half monochromator width) /i? = 
±250/2000 = ±12.5%. This is acceptable because we do 
not aim for fluU-flcdged structure determination but only 
for monitoring structural changes.-^ 

G. Instrument Housing with Argon Filling 

The instrument housing has a surface of about 36 m^, 
and a volume of nearly 60 m'^. Its walls and ceiling are 
made of steel plates of 5 mm thickness, anchored to a steel 
frame, and proven to be gas tight. Towards the inside, a 
100 mm polyethylen layer moderates fast neutrons. The 
innermost cladding of walls and ceiling consists of Cd 
plates of 1 mm thickness. 

The housing is reasonably gas tight and equipped 
with pneumatics and controls so that it can be filled 
with argon. This reduces losses by air scattering in the 



secondary spectrometer, where neutrons have to travel 
about 4.3 m through atmosphere. In the primary spec- 
trometer, another 0.3 m are in atmosphere; the remaining 
fligh path is in a vacuum chamber that prolongates the 
chopper casing towards the monochromator, as indicated 
in Fig. [H Displacing 90% of the air by argon leads to an 
intensity gain of more than 35%, in good accord with 
expectations from tabulated neutron cross sections. For 
technical details of the argon filling, see the Supplemental 
Material.-^ 



H. Sample Environments 

The standard sample environment, used for most ex- 
periments, is a cryostat SIII-950T (Janis) It is a closed- 
cycle refrigerator system with cold head and compres- 
sor manufactured by Sumitomo Heavy Industries. In the 
low-temperature range from 4 K to 320 K, the sample is 
cooled via exchange gas. In the high-temperature range, 
from 290 K to 650 K (700 K for limited amount of time), 
the system acts as a cryofurnace, with the sample in vac- 
uum. The sample tube has a diameter of 60 mm; samples 
are allowed to have diameters of up to 55 mm. 

The standard sample stick is equipped with Si diodes 
DT670A for temperatures from 4 K to 500 K. A high- 
temperature stick with Pt 1000 thermometers covers 
temperatures from 10 K to 700 K. Gas- loading equip- 
ment is under preparation. 

Sample environments are inserted and operated from 
the platform on top of the instrument housing. If the 
tail diameter exceeds 160 mm, the secondary spectro- 
meter must be temporarily realigned to a sample po- 
sition slightly out of the focus of the monochromator. 
To change the sample environment, it is unavoidable to 
break the argon atmosphere. Therefore, experiments 
with non-standard sample environments are preferen- 
tially scheduled at the beginning or at the end of a reactor 
cycle. So far, experiments have been performed with a 
dilution cryostat that gives access to temperatures down 
to 30 mK. Other environments will be adapted and tried 
as requested by users. 



IV. INSTRUMENT CONTROL AND DATA 
ACQUISITION 

A. High-Level User Interface 

During regular experiments, physical access to the in- 
strument is only needed for changing samples. All other 
operations can be performed remotely. A graphical user 
interface (GUI) allows 

• to open and close the instrument shutter, 

• to start and stop the monochromator drive, and to 
set its velocity amplitude. 
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• to set the sample temperature, and to start tem- 
perature ramps, 

• to start and stop the data acquisition, to save files, 
and to reset histograms, 

• to submit, edit, start and stop experiment scripts. 
It displays 

• warnings and error messages, 

• the status of data acquisition, sample environment, 
and peripheral systems, 

• current count rates for all detectors, separately for 
direct, elastic, and inelastic scattering, 

• the latest acquired spectrum, 

• the latest acquired diffractogram, 

• the time evolution of various parameters in the last 
hours or days, 

• the status of experiment-script execution and the 
script commands ahead. 

As the instrument is operated from different work- 
places (cabin in the neutron guide hall, office, home), 
there is no simple way to maintain a traditional log book. 
Therefore, we have replaced paper by digital storage, us- 
ing a simple wiki system~2, as our instrument log. Besides, 
there are automatically generated event and parameter 
logs. 

B. Instrument-Control Daemons 

As a leading principle in the design of our instrument- 
control software, we requested that no neutrons should 
be lost because of peripheral software or network prob- 
lems. Therefore, the instrument's subsystems are con- 
trolled by a number of independent background processes 
(daemons) that continue to work even if the connection 
to upper-level processes is interrupted. 

The following daemons are the most important ones: 

• spectral data acquisition (Sect. ITVDt . 

• diffraction data acquisition (Supplemental 
Material >2S 

• slow control daemon (Sect. IIV E[) . 

• temperature-controller controller, supervising the 
sample-environment temperature controller, 

• Doppler driver, supervising the control PC supplied 
by the manufacturer, 

• event logger, receiving log lines from other daemons 
and writing them to an event log, thereby docu- 
menting user interventions, experiment-script exe- 
cution steps, system warnings and error messages. 




IRQ 
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FIG. 5: The data acquisition VME crate.— The signal pro- 
cessor on the SIS 3100 sequencer board is a ADSP21061. An 
optical fiber connects the crate to a SIS 1100 PCI card on the 
dedicated data-acquisition PC. 

• parameter logger, querying every 30 s some fifty 
status parameters from other daemons and writing 
them to a parameter log, 

• documentation server and raw-data postprocessor 
(Sect. lVA| . 

• script engine, executing user-written experiment 
scripts, 

• GUI server, through which data and commands are 
channeled to and from the graphical user interface. 

The spectral acquisition is time critical; therefore it is 
written in the programming language C, and it is running 
on dedicated computer. For all other processes, computa- 
tional speed is less important than clarity of code; there- 
fore they are written in an expressive scripting language. 
Ruby. 

All computers are connected by ethernet (instrument 
LAN). Interprocess communication (IPC) is based on 
TCP (Transmission Control Protocol) messages. Each 
daemon is listening on a TCP port for commands given 
by other daemons. For instance, the spectral acquisition 
daemon receives commands (start, stop, save, reset) from 
the GUI server, and status queries from the parameter 
logger. 

C. Spectral Acquisition Electronics 

The data acquisition is done in a VME crate (Fig. [S]) 
using components and software from Struck Innovative 
Systeme (SIS).'^^ The input consists of TTL pulses from 
the neutron-counter preamplifiers and from the optical 
encoders sensing the translation of the monochromator 



7 



Doppler drive (200 ticks per mm) and the rotation of the 
PST chopper (20480 ticks per revolution). The Doppler 
and chopper ticks are counted up and down in home- 
made counter boards, then they are transmitted to a 
32 bit latch. The neutron counts are transmitted via 
a trivial home-made interface board to a 32x32 scaler. 

The sampling period is r = 20 ps. Once every r, a dig- 
ital signal processor reads the data from the latch and 
from the scaler, writes them to DRAM, and resets the 
scaler. The data form a fixed-length block of 20 32-bit 
words. In the following, each block shall be identified by 
an index i, corresponding to a time ti — const -I- it. One 
word encodes the Doppler position Zi and the chopper 
phase Xi (16 bits each); the remaining 19 words riji ac- 
comodate the neutron counts from up to 19 detectors j. 
The DRAM is divided into two buffers with I'OOO'OOO 
words each. It takes 50'000 data blocks or 1 s to fill one 
buffer. While one buffer is being filled, the other buffer 
is read out via an optical link by the upstream PC. 

D. Spectral Acquisition Daemon 

The data-acquisition PC (Opteron 2x2.4 GHz, 64bit), 
is dedicated to running the spectral acquisition daemon. 
It accumulates and stores histograms containing neutron 
counts per energy channel and per detector. These raw 
spectra are only a few correction steps (Sect. IV Bp away 
from the physical scattering law S{q,uj). 

The daemon consists of an endless loop that waits for 
one of the two DRAM buffer on the VME crate to be 
filled. As soon as this is the case, its contents is read 
via the optical link and pushed on a stack. When the 
stack covers at least one Doppler oscillation period, it is 
processed. A sine function with four free parameters (off- 
set, amplitude, frequency, phase) is fitted to the Zi, using 
starting values retained from the previous fit. Then, a 
loop over the data blocks i is performed, to be described 
in the next three paragraphs. 

Depending on the chopper phase Xiy the data block is 
assigned a category c(xi), which can be indirect, direct, 
or ambiguous, depending on whether or not the neutrons 
have travelled the detour sample-analyzer-sample or not. 
In Sect. IVlO we will show how this categorization is set 
up. Ambiguous data blocks are not processed further. 

To assign our scattering events to an energy channel lo, 
we need to know the time 

U = - (wf + v{U)){L° - z{U)) - viL'i (4) 

at which the neutron has been backscattered by the 
monochromator. Here, ti is the detection time, and Li^f 
is the path travelled before and after scattering by the 
sample. The path Lf depends on the category c: it is 
only about 10 cm in the direct case, but 4 m more in 
the indirect case. The path offset z{ti) and the velocity 
modulation v(ti) are given by position and speed of the 
monochromator at time ti. This makes Eq. Q nonlin- 
ear; it is most easily solved by iteration. For a given ti. 



the analytic derivative of the fit function z(t) is used to 
obtain v{ii). Eq. ^ then yields uj{vi). A quick estimate 
shows that the delay ii — ti is not negligible: Neutrons 
in the indirect category travel more than 6 m from the 
monochromator to the detector, which takes about 10 ms, 
corresponding to up to 1/10 of a Doppler period. 

On starting the data acquisition for a given Doppler 
velocity amplitude, an equidistant uj mesh has been 
chosen, with channels [ujk — Aw/2,ti;fc + Auj/2]. This 
makes it easy to determine the pertinent channel k{i) = 
[{uj{vi) — {uq — Aw/2))/AwJ. Finally, the neutron his- 
tograms of the pertinent chopper category c{xi) are in- 
cremented: Hji^ Hjj. + nji for all detectors j. The time 
spent in energy channel k(i) is incremented accordingly: 
Nj: <— Nj: + 1. Here ends the loop over i. 

All these computations must be completed before the 
next DRAM buffer is overwritten. After some optimiza- 
tions, this requirement could be fulfilled with sufficient 
margin: In years of operation, not a single data block has 
been lost because of computer overload. As a side prod- 
uct, a self-contained C library for Levenberg-Marquardt 
least-squares minimization has been released as an open- 
source project^ serving practitioners from all areas of 
science and engineering. 

Our histogramming relies on the assumption that the 
motion of the monochromator accurately follows a sim- 
ple analytic function. Analyses of recorded Zi{ti) traces 
have shown that the motion accomplished by the Aero- 
las Doppler drive is indeed almost perfectly sinusoidal, 
deviations being negligible on the scale set by the energy 
resolution of the spectrometer. 

The accumulated histograms are written to files when 
a preset time interval has elapsed or when a save com- 
mand is received through the daemon's IPC interface. 
Files are saved to the local disk, which is remote read- 
able through NFS (Network File System). The raw- 
spectra, as other data produced by SPHERES, are for- 
matted as self-documenting YAML files. YAML is a 
data-serialization specification with intrinsic support for 
basic data structures like arrays and hashes and with an 
emphasis on human readability; parsers are available for 
all major programming languages.— 



E. Safety Instrumentation and Slow Control 

SPHERES is equipped with a fail-safe programmable 
logic controller (Fehlersichere speicherprogrammierbare 
Steuerung F-SPS, Siemens Simatic S7 300F) for safety 
instrumentation and slow control4i. The F-SPS is con- 
nected via internal bus or Profibus-DP to various sensors 
and actors. A touchpanel next to the instrument door 
provides a status display; it provides the only regular 
user interface for safety-critical subsystems like the door 
lock and the argon pneumatics. 

A simple daemon running on a compact-PCI computer 
couples the F-SPS to the instrument LAN. It has read 
access to all sensors and write access to safety-uncritical 
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actors, including the shutter opener, which can be acti- 
vated if the instrument housing is closed and locked and 
a manual switch is in "remote control" position. 

V. DATA PROCESSING 
A. Raw-Data Postprocessor 

The raw files produced by the spectral data-acquisition 
daemon contain the energy mesh w^, the histograms Hjf, 
and N^, some information about the instrument and 
data-acquisition setup, and the date and time when the 
measurement was started and stopped. Had all other 
daemons crashed, this would be sufficient to allow a valid 
data analysis. In regular operation, however, users ex- 
pect data files to contain additional information like ex- 
periment title, sample name, and sample temperature. 
Since the spectral acquisition daemon is deliberately kept 
ignorant about these things (not essential for its func- 
tioning, and potential sources of bugs and failures), the 
wanted information must be added ex post. 

This is done by the raw-data postprocessor. Running 
on another computer than the spectral acquisition dae- 
mon, it reads a raw-spectrum file via NFS, and looks up 
the start and stop time entries. Then, it opens the pa- 
rameter log, reads the lines recorded during the spectral 
measurement, and computes mean value and standard 
deviation of the sample temperature and of some other 
parameters. Textual information like experiment title 
and sample name, usually entered by the user via the 
GUI, is retrieved from another dedicated file system. All 
this is written to a user-friendly aggregated raw-data file. 

B. Raw-Data Reduction with SLAW 

While some users want to work on the raw data, most 
ask for standard software to deliver normalized scattering 
laws S{q,u!). To provide this service, a new raw-data re- 
duction program Slaw^ has been developed. Compared 
to legacy software like SQW of the ILL, Slaw offers bet- 
ter support for four-dimensional data sets S{q,uj;T) or 
S{q,io;t) obtained in time-resolved inelastic scans. Also, 
Slaw does not bind users to a specific data format; it 
supports a variety of output formats, and new output 
routines can be easily added ad hoc. 

Slaw is controlled by an input script that contains 
configuration commands and a map that assigns raw 
scans to output files. Slaw processes the raw data 
through the following steps of which all but the first and 
the last are optional: 

• read raw data files, adding or subtracting subscans 
as needed; 

• perform binning in energy, scattering angle, or time 
slots, possibly eliminating some of these dimen- 
sions; 



• subtract empty-can measurement, optionally 
weighted with a self-absorption factorji^ 

• normalize neutron counts Hj/^ to the measuring 
time Nj:T spent in a given energy channel; 

• correct for incident flux variations (Sect. IV C|) : 

• divide counts by elastic intensity of a reference mea- 
surement (usually vanadium or a low-temperature 
measurement with the current sample); 

• save the scattering law. 

Slaw has alternative input routines to process raw 
data from other backscattering spectrometers. Its struc- 
ture is generic enough to allow future extension to time- 
of-flight spectrometers. 



C. Instead of a Monitor 

Most inelastic neutron spectrometers have a high- 
transmission neutron counter, located a little upstream of 
the sample, to monitor the incoming flux. In SPHERES, 
the windows of such a monitor would cause undesirable 
background. We therefore prefer an indirect normaliza- 
tion. 

We must account for two kinds of fiux variation: tem- 
poral variations of the reactor power, of the filling of the 
cold source, and of the monochromator orientations of 
the upstream instruments, and spectral variations of the 
primary spectrometer transmission. The temporal vari- 
ations of the incident beam can be measured with good 
accuracy without inserting any device into the beam: it is 
sufficient to monitor the 7 radiation emitted by the alu- 
minum windows at the end of the convergent guide and 
at the chopper entry. Such a 7 monitor, still missing, will 
be installed soon. 

The spectral variation of the flux at the sample is de- 
termined by spectrum of the cold source, by the neu- 
tron guide system, by the velocity selector, and most 
importantly, by the PST chopper and the monochroma- 
tor. To make things more complicated, the angle of the 
Bragg deflection in the PST chopper is energy depen- 
dent, and the monochromator is moving in space so that 
the geometric-optical focussing onto the sample is imper- 
fect, with a non-trivial interplay of energy, direction, and 
lateral spread. To account for all this, the spectral dis- 
tribution of the incident flux is best measured either by 
placing a neutron counter exactly at the sample position, 
or by counting neutrons that are scattered from the sam- 
ple directly into the regular detectors. The latter method 
is regularly used as the "pseudomonitor" of SPHERES. 
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VI. PERFORMANCE supermirror coating of the focussing neutron guide is fol- 

lowed by a conversion to fast neutrons in the guide's lead 
A. Flux and Count Rates shielding. 



The incident flux at the end of the convergent neutron 
guide has been measured at two occasions with gold foil 
activation. The results were (1.85 ± 0.05) • 10^ cm~^s~^ 
in November 2005, and (1.19 ± 0.05) • 10^ cm'^s"! in 
February 2012. A decrease of the order of 10% was to be 
expected from the two breaks in the upstream neutron 
guide fSect. UlL^ . The remaining loss is not understood; 
it must be feared that the neutron guides have degraded 
substantially. To detect further degradation and to deter- 
mine necessary action, it will be important to remeasure 
the flux more regularly in the coming years. 

The flux at the sample position is more difficult to mea- 
sure; for details see the Supplemental Material.^- The 
best current estimate is ^sam = 1.8 • 10^ s~^ within a a 
cross-section of 25 x 40 mm^. Here and in the follow- 
ing, fluxes are reported for the open chopper state. The 
average flux seen by a gold foil is only half of it. 

For typical samples made for transmissions between 
80% and 90%, elastic count rates (/)ei are of the order 
1000 s~^. Here as in the instrument's GUI display, such 
numbers are per detector, and under the assumption that 
all time is spent counting into this single channel. The 
pseudomonitor rate 0dir of neutrons that are directly 
scattered from the sample into the detector, is about 
twice as large. Reasons why 4'e\/4'di-[ is smaller than 1 
include: solid angle mismatch (the detector block has 
a larger solid angle than the analyzers, to ensure that 
all backscattered neutrons have a chance to reach a de- 
tector), absorption and scattering losses in the sample 
(which must traversed a second time by backscattered 
neutrons), losses along the 4 m from the detector to 
the analyzers and back, imperfect energy selection in the 
monochromator, finite backscattering probability in the 
analyzers. 



B. Dark Counts 

Dark count rates in each of the large-angle detectors 
are of the order of 0.02 s~^ when the reactor is off, and 
0.1 s~^ when the reactor is at 20 MW and the intru- 
ment shutter closed. This is negligible compared to the 
background produced in our primary spectrometer. 

When the instrument shutter is open and the chopper 
is running, but the exit window of the chopper (just be- 
fore the sample position) is closed by a Cd mask, then the 
dark count rates in the large-angle detectors are about 
0.53 s^^. Only a small part of this background is due 
to slow neutrons that penetrate from the primary into 
the secondary spectrometer, circumventing the chopper 
exit window. Most of the background is independent 
of the chopper phase, which indicates an origin between 
the shutter and the chopper. This has been confirmed 
by experimental tests. Probably, a n-7 conversion in the 



C. Chopper Modulation of Count Rates 

Along with the duty-cycle chopper, a backscattering 
spectrometer needs a tinie-of-flight logic to discriminate 
neutrons that are backscattered by the analyzers from 
neutrons that are scattered directly from the sample into 
the detectors. In SPHERES, this time-of-flight discrim- 
ination is performed by the spectral acquisition daemon 
fSect. irVDj) . based on the chopper phase reading that is 
saved every 20 \\s along with the Doppler position read- 
ing and the neutron counts. 
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FIG. 6: Neutron counts versus chopper phase. Data from the 
same low-temperature resolution measurement as in Fig. [7] In 
the direct phase, neutrons are scattered without energy anal- 
ysis from the sample into the detectors. In the indirect phase, 
registered neutrons come from the analyzer. If the scattering 
is purely elastic, as in this example, then the inelastic signal 
is just undesired background. 

The parametrization of the discrimination must be 
determined empirically. A chopper-phase histogram is 
saved in each raw data file. This histogram is based 
on a rectangular grid with 120 chopper-phase channels, 
but only two energy channels, to distinguish roughly be- 
tween elastic and inelastic scattering. Fig. |6] shows such 
data from a resolution measurement. There is a range 
of nearly 60° with very high count rates, independently 
of the Doppler velocity, due to direct scattering. There 
is another range, labelled indirect, where almost all neu- 
trons have been backscattered by an analyzer, as can 
be inferred from the strong difference between the elas- 
tic and the inelastic channel. And there are intermediate, 
ambiguous ranges where directly scattered neutrons form 
an acceptable background in the inelastic channels. Part 
of this transitory range is due to the finite time it takes 
the chopper to open and to close. Other contributions 
include the following: 

• In the present chopper, the graphite crystals are 



10 



kept in place by aluminum corner pieces. When 
these pieces cross the primary neutron beam, they 
cause incoherent scattering into 47r, thereby send- 
ing a diffuse neutron cloud into the secondary spec- 
trometer. 

• Some directly scattered neutrons are retarded 
by phonon downscattering, causing a sample- 
dependent extension of the chopper-closing phase. 

• Multiple scattering (first scattering in the sam- 
ple, second scattering by the sample environment 
or somewhere in the secondary spectrometer) also 
presents a sample- dependent contribution to the 
extended chopper-closing phase. 

To parametrize the data acquisition daemon, the in- 
strument responsible must make a relatively arbitrary 
choice which range of chopper phases shall be considered 
as indirect. In general, users of SPHERES are more inter- 
ested in an excellent signal-to-noise ration than in maxi- 
mum count rates. Therefore, a relatively strict chopper- 
phase discrimination that favors signal-to-noise ratio on 
the expense of total count rates is preferred. For most 
measurements in the past years, we have set the indirect 
interval to 35°, which means that the duty cycle (the 
fraction of time actually used for incrementing scatter- 
ing histograms) is no more than 35/120 = 0.29. 

For the future, one might consider a small modification 
of the data acquisition code to simultaneously record en- 
ergy histograms for different chopper-phase discrimina- 
tion settings. Users could then decide ex post whether to 
analyze a "strict" data set with optimum signal-to-noise 
ratio, or a "loose" data set with higher count rates. 



D. Resolution and Signal-to-Noise Ratio 

The instrumental resolution can be empirically deter- 
mined by measuring the elastic scattering of a solid. The 
scattering should be predominantly incoherent to guaran- 
tee sufficient intensity at all scattering angles. Vanadium, 
the elastic incoherent reference scattering par excellence, 
is not well suited because its absorption cross section at 
A = 6.27 Ais more than three times larger than its scat- 
tering cross section. Therefore, we prefer a simple plastic 
foil (Hostaphan RN 250: a polyethylene terephthalate foil 
of 0.25 mm thickness) as a our reference scatterer. 

In most experiments, however, the resolution is mea- 
sured on the user sample, cooled to a low temperature 
where almost all scattering is elastic. This has the ad- 
vantage that the sample geometry is exactly the same as 
in the production measurements. 

Fig. [7] shows such a resolution spectrum from a user 
experiment. The Doppler drive was running with the 
maximum velocity amplitude of 4.7 m/s, resulting in a 
dynamic range of ±31 |J.eV. To determine the resolution 
width and the elastic amplitude, the centre of the elas- 
tic peak (|fta;| < 0.6 M-eV) has been fitted by a Gaussian, 
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FIG. 7: A resolution measurement with a user-provided sam- 
ple (a powder of hydrogen-rich organic crystals at 3 K, in 
a flat cell calculated for 90% transmission). Cumulated data 
from the six large-angle detectors with the best signal-to-noise 
ratio. The solid line is a Gaussian fit to the central region 
< 0.6 ^eV. 



which is an excellent approximation though there are pro- 
nounced non-Gaussian wings at larger energy transfers. 
The full width at half maximum (fwhm) is obtained by 
multiplying the standard deviation of the Gaussian with 

When the Doppler drive is run with small to moder- 
ate velocity amplitudes, the resolution at large scatter- 
ing angles is typically between 0.62 and 0.65 |J.eV. At 
the maximum Doppler frequency used here, the resolu- 
tion is slightly reduced (to 0.66 M-eV), probably because 
of distortions of the monochromator support. 

The signal-to-noise ratio (SNR) depends obviously on 
how we define the "noise" . If we average only over en- 
ergies far off the central peak, say \hu!\ > 20 |J.eV, then 
we obtain in the present example 1750:1. The SNR also 
depends on the sample: if the sample is too thin, then 
the sample-independent constant background contributes 
relatively more; if the sample is too thick, then many 
backscattercd neutrons are lost when they retraverse the 
sample. 



VII. COMPARATIVE STATUS AND 
PERSPECTIVE 

In Table |lTl some design and performance parameters 
of the three backscattering spectrometers IN16, HFBS 
and SPHERES are compared. These parameters are in- 
tercorrelated as follows. 

The instrumental resolution is within 5% proportional 
to the chosen thickness of the Si(lll) monochroma- 
tor /analyzer wafers. This confirms the fundamental im- 
portance of the crystal choice for the instrument per- 
formance, and it shows that all three instruments very 
closely attain the optimum resolution allowed by their 
crystal sets. 
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instrument 


IN16" 


HFBS 


SPHERES 


reactor power (MW) 


58 


20 


20 


PG(002) speed (m/s) 


75 


250 


100 


!3i(lll) tluckness (vim) 


700 


750 


500 


resolution fwhm ((ieV) 
flux at sample'' (10*^ s"^) 


0.8-0.9 


0.93 


0.62-0.65 


1.2 


2.4 


1.8 


rescaled flux'^ 


0.17 


0.91 


1 


signal-to-noise ratio'' 


1000' 


600 


1750 



"With the unpolished, deformed monochromator/analyzcr set 
that is used in almost all experiments. 

''While chopper open. Au foil activation gives half this value. 

■^Flux at sample, divided by resolution width and reactor power, 
arbitrarily normalized to the value of SPHERES 

''Best value reported for a non-absorbing sample of optimum 
thickness. 

■^Not mentioned in the instrument descriptionsi^^ii^i^, but con- 
firmed by instrument users. 

TABLE II: Comparison of the three backscattering spec- 
trometers IN16 of the ILL^-'i^ HFBS at NIST.^i^ and 
SPHERES of the JCNS at FRM II. 

Ceteris paribus, the flux at the sample should be pro- 
portional to the reactor power. In a very good approx- 
imation, it should also be proportional to the width of 
the energy band selected by the monochromator. There- 
fore, the figure of merit that describes the overall effi- 
ciency of the cold source, of the neutron guides, and 
of the primary spectrometer is flux divided by reactor 
power divided by resolution width. The comparison of 
this rescaled flux shows how much the compact design of 
HFBS and SPHERES, due to Schelten and Alefeld," is 
superior to the double-deflector layout of IN16. 

The primary spectrometer transmission is also some- 
what related to the signal-to-noise ratio: to maximize the 
latter, we installed in SPHERES ad hoc several slits that 
cut away some 10% of the incoming and of the backscat- 
tered beam in order to prevent neutrons from being trans- 
mitted through the closed chopper. The next chopper, 
currently under construction, will hopefully allow us to 
take away these slits and to increase the flux at the sam- 
ple accordingly. Furthermore, in this new chopper the 
PG(002) deflector crystals will be moved with about op- 
timum speed, which will result in yet another important 
gain in flux. 

At small scattering angles, analyzers cannot be aligned 
to exact backscattering geometry because detectors must 



be placed outside the incoming neutron beam. This 
makes the resolution much wider and more assymetric 
than in the regular backscattering detectors. Supported 
by simulations, we have improved the resolution at small 
angles by breaking the azimuthal symmetry of the ana- 
lyzer rings, as will be described separately. 
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